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An investigation of possibilities for superconductivity me- 
diated by spin fluctuations in some elementary metals is mo- 
tivated by the recent discovery of superconductivity in the 
hep high-pressure phase of iron. The electronic structure, 
the electron-phonon coupling (X p h) and the coupling due to 
spin-fluctuations (X s f) are calculated for different phases and 
different volumes for four elementary metals. The results 
show that such possibilities are best for systems near, but 
on the non-magnetic side of, a magnetic instability. Fee Ni, 
which show stable magnetism over a wide pressure range, is 
not interesting in this respect. Ferro- and antiferro-magnetic 
fluctuations in hep Fe contribute to a relatively strong cou- 
pling in the pressure range where superconductivity is ob- 
served. The absence of fluctuations at large q-vectors makes 
fee Pd only moderately interesting despite its large exchange 
enhancement for q=0. Fee Co at high pressure (~ 0.5 Mbar) 
behaves as an improved version of Pd, where the fluctuations 
extend to larger q. The estimations of To, which reproduce 
the experimental situation in Fe quite well, suggest a mea- 
surable To for the high-pressure phase of fee Co, while the 
estimate is lower for the ambient-pressure phase of fee Pd. 
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I. INTRODUCTION. 

Superconductivity has recentW been reported in the 
hep high-pressure phase of ironEru. The maximum of the 
superconducting transition temperature Tc, about 2 K, 
is in the middle of a pressure range of 150-300 kBar where 
superconductivity has been observed. Superconductivity 
is traditionally explained from electron-phonon coupling. 
This mechanism is sensitive to magnetic interactions and 
it would not be expected to find superconductivity in 
magnetic or almost magnetic materials. Alternatively, 
superconducting pairing can be mediated by spin fluctua- 
tions, and it might be that hep Fe should be explained by 
this type of mechanism, since anti-ferromagnetic (AFM) 
instabilities are likely in the interesting pressure region. 
Further, some magnetic or nearly magnetic conrnonnds 
have also been found to be be superconductingem, al- 
though the Tc's are not very high. Both the electron- 
phonon coupling and ferromagnetic (FM) instabilities re- 
quire a large electronic density-of-states, and the question 



is if superconductivity in such cases is the standard one 
caused electron-phonon coupling where Tc is limited by 
spin fluctuations, or if it is mediated by spin fluctuations 
in equal spin pairing (ESP). The fact that the supercon- 
ductivity is found only in very pure samples is an indi- 
cation that mechanism is the latter one, since impurity 
scattering is expected to interfere with ESP in a similar 
manner as magnetic impurities interfere with pairs which 
are hold together by electron-phonon coupling. In addi- 
tion, theory indicates the pressure dependence of Tc in 
Fe is inconsistent with the mechanism of electron-phonon 
couplingB'cl. The recent discoveries of magnetic supercon- 
ductors led to a renewed interest in theories of calculating 
the coupling strength, A s / , for spin fluctuations ancLhow 
such fluctuations can mediate superconductivityoEj. 

Superconductivity in some borocarbides, containing 
magnetic xare-earth elements is also surprising and 
interestingcj, but the situation is probably quite different 
from Fe and the compounds mentioned above, since the 
magnetic f-electron of the rare-earth atom is very local- 
ized and electronically isolated from the itinerant valence 
electrons. High-Tc copper oxides paght be superconduc- 
tors because of AFM nuctuationsli3, but stripe-like spin 
waves may interact with phononstj, and the mechanism 
not yet clear. 

This work consider itinerant electrons in pure metals 
where magnetic fluctuations and superconductivity must 
coexist. Near-magnetic pure elements are particularly 
interesting in studies of superconductivity and magnetic 
fluctuations, both experimentally and theoretically. Non- 
compound elements can be made very pure without in- 
trinsic defects caused by non-stoichiometry or vacancies 
in sub-lattices of compounds. This allows for ESP to de- 
velop without much interference from impurity scatter- 
ing. Theoretically it is also preferable to consider pure 
metals over compounds. The electronic structure is sim- 
pler and especially for studies of extended AFM fluctu- 
ations by the method of frozen spin waves, it avoids the 
complication of having several types of sites where the 
local magnetic response can be different. In this work 
we study four metals, Fe, Co, Ni and Pd, which all are 
known to be ferromagnets or strongly enhanced param- 
agnets. Different phases are considered, and total energy 
calculations are used to determine the stability of hep, 
fee and bee phases as function of pressure. A previous 
report considered FM and AFM spin fluctuations sepa- 
rately in hep Fe, and it was concluded that the combined 
effect could explain the observed Tea This study is here 
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extended to include the other three metals and where a 
weighting of different fluctuations are made, similar to 
what is made for phonons. 

II. THEORY. 

The electronic structures are calculated using the lin- 
ear muffin-tin orbital method in the local spin density 
approximation (LSDA)Elc£l and the generalized gradi- 
ent approximation (GGAo) scheme, which includes non- 
local corrections to the potential. The basis set includes 
states up to £=3. The c/a-ratio in the hep structures 
is fixed to f.59 for all volumes. This is close to the 
calculated value within a wide pressure range for FeEa, 
and a few tests with other ratios show that this choice 
is good also for Co. The number of k-points vary with 
the structure. For fee and bec it 505 and 506 points in 
the irreducible part of the Brillouin Zone (BZ), for hep 
252, and for the supercells with long wavelength AFM 
configurations we use the same density of points. The 
number of points is close to 500/n where n is the num- 
ber of atoms per cell. The use of GGA instead of LSDA 
potential is known to be crucial for a correct description 
of the fee vs. bec stability of Fe, but other properties 
are not very different in the two types of density func- 
tional potentialsEI A previous report on Fe used mainly 
LSDA for the calculations of A s /H Here, these results are 
complemented with calculations using GGA for the most 
interesting configurations in hep iron. The magnetic cou- 
pling constants are similar (within-^ 5 percent) to those 
calculated in the LSDA. As in ref.EiJ it is found that two 
AFM configurations compete with the NM state in hep 
iron at some pressures. The first (AFM-I) has opposite 
polarization on each alternate z-plane layers, while the 
second (AFM-II) has opposite polarization on each layer 
perpendicular to the x-axis. The results for Ni, Co and 
Pd are new, and only results from GGA calculations are 
reported. 

The electron-phonon coupling parameter X ep is calcu- 
lated from the electronic structure; 

A ep = NV 2 /Mu 2 = NV 2 /K (f ) 

where N is the density of states at Ep, M the atomic 
mass, u a phonon frequency and K a force constant. The 
latter is related to the change of total energy E due to 
an atomic displacement r, K = d 2 E/dr 2 . The Hopfield 
parameter NV 2 contains V, the Fermi surface (FS) aver- 
age of the matrix element < ipk{dv(r) / 'dr)ipk >, the first 
order change in energy due to the displacement, where 
ipk is the wave function at the k-point k, and v(r) is the 
potential. For elementary metals it is a good approxima- 
tion to calculate, V within the rigid muffin-tin approxi- 
mation RMTAB, which in the case of LMTO band re- 
sults is slightly— modified to the rigid Wigner-Seitz (WS) 
approximation^. The alternative to calculate V from 
the change in band energies, when the structure is 



deformed as in a 'frozen' phonon, demands a larger com- 
putational effort since several large supercell calculations 
are needed to map out representative phonons. The force 
constant is fitted to the experimental value of the De- 
bye temperature Od at zero pressure, and then .scaled to 
other volumes by the calculated bulk modulusEj. This 
scaling is assumed to be valid also between the hep and 
bec phases for Fe and between the hep and fee phases of 
Co. The calculated values of A ep at equilibrium volumes 
agree well with independent results for the elements con- 
sidered hereESnJ. This is despite the use of Wigner-Seitz 
cells instead of non-overlaping Muffin-Tin spheres in the 
evaluation of the matrix elements. 

Interesting spin fluctuations occur near magnetic in- 
stabilities when only a small energy is needed to excite a 
magnetic wave in the material. The energy is zero at the 
magnetic instability and it can be very low close to the 
instability, lower than the energy of typical phonon en- 
ergies. It is also noted that spin fluctuations can be the 
result of a divergence of the generalized susceptibility, 
x(q,co), implying that Fermi surface nesting for critical 
values of q is partly responsible for the divergence. This 
is also the case for phonon softening when some g-modes 
become unstable. Phonons with q outside the BZ have 
no sense, since the wave length cannot be shorter than 
two interatomic distances. Spin- or charge- density waves 
with wavelengths shorter than the size of an atom can be 
imagined, but if they were important they should prob- 
ably exist already on the level of a single atom, before 
the scattering within several sites in a lattice can set up 
a FS and nesting. Since the structure always seems to 
be important for the properties of a material, we assume 
that fluctuations corresponding to g-vectors outside the 
BZ are not important. Thus, we consider only "soft" 
spin fluctuations with energies and wavelengths which 
are similar as for phonons. If adiabatic conditions and 
the Migdal theorem are valid for phonons they should 
also be valid for soft spin fluctuations. _ 

As was pointed out in a previous workQ, the analogies 
between soft spin-wave excitations and phonons can be 
utilized. A spin wave with a moment m and a phonon 
with lattice-displacements "r" both generate perturba- 
tions of the potential, although the perturbations are 
of different kinds. Magnetic moments and the spin- 
polarized perturbation of the potential have spherical 
symmetry within each site, except for very short wave 
lengths of the spin waves. The matrix elements V are 
mainly monopolar, which couple electronic states of equal 
£. Atomic displacements generate mainly an asymmet- 
ric (and non-polarized) perturbation of the potential at 
each site. Ionic compounds, where the displacements are 
associated with large, Madelung shifts, can have some 
monopolar couplingEZI. But here, in neutral pure ele- 
ments, the matrix element V is dominantly dipolar with 
the selection rule I — £' = ±1. In addition there is a tech- 
nical difference of how the frozen phonon and spin wave 
calculations are made. By applying a magnetic field £ 
one obtains a magnetic moment m in frozen spin wave 
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configurations. However, the displacements in frozen 
phonon calculations are imposed instead of being a re- 
sult of an applied force. This difference would disappear 
if we did so-called "fixed-spin- moment" (FSM) calcula- 
tions instead of applying the magnetic field. The total 
energy in such case is the free energy for the configura- 
tion with the fixed moment, while with imposed field the 
free energy is the total energy plus the term field times 
moment. In practical calculations we prefer the method 
with applied fields which is simple to use both for FM 
and AFM cases. 

The calculation of mass enhancement due to FM or 
AFM spin fluctuations, A s /, is analogous to that for 
phonons, but with "r" replaced by "m". 

\ sf = NV 2 /K = N < ijj k (dv(r)/dm)ij k > 2 /(d 2 E/dm 2 ) 

(2) 

Calculations for spin wave excitations with general q- 
vectors require supercells, as in the case of frozen phonon 
calculations. The calculations are made for two differ- 
ent configurations, one with an applied field and one 
without, to give two band structures and two magne- 
tizations. If the changes of the free energy (E) are har- 
monic, i.e. quadratic functions of the induced magnetic 
moment, we have Eg — Eq + Jm 2 where m is the mo- 
ment (induced by the field £) per atom. The constant 
J = \d 2 E/dm 2 = (Eg - E Q )/m 2 is calculated from spin 
polarized results, where the applied fields range from 1 to 
10 mRy. The local Stoner enhancement, S, is defined as 
(eg — eo)/£, where eg is the increase of the exchange split- 
ting at an atom (obtained from the logarithmic deriva- 
tive of the d-band) induced by the field £. (K-dependcnt 
Stoner enhancements can be obtained by taking the ex- 
change splitting from the band energies at each k-point.) 
The S- values are generally smaller for large £, indicating 
some aversion against large moments and non-harmonic 
variations of J as function of field. 

The monopolar matrix element < ipkdvtpk' > can be 
calculated directly from the change in the spin part of 
the potentials between the two cases, with and without 
field. This permits calculation of the anisotropy of the 
coupling, X(k, k'), which is needed in Eliashberg theory 
to determine if the superconductivity is s-, p-wave or of 
some higher symmetry. In this work we do not determine 
the anisotropy, but only the FS average of X(k,k). In 
this case it is convenient to determine the matrix element 
directly from the change in band energies, simultaneously 
to the 'frozen spin wave' determination of the change in 
total energy. Hence, < ipkdmpk > 2 = (e™ — e") 2 , where 
is the band energy at point k for the configuration 
belonging to m. In the harmonic limit we obtain finally; 

X sf =N< (el! - el) 2 > FS /2(E m - E ) (3) 

Here < >fs means a FS average. The matrix element is 
more efficient for FM than for AFM fluctuations, because 
a FM field applied in a system where one band is domi- 
nant at Ep, will split all bands about equally. The same 



difference (e™ — e°) will appear almost everywhere over 
the FS. For induced AFM configurations it is more diffi- 
cult to visualize these differences. Positive and negative 
differences are likely, which means that small differences 
should exist at some sections of the FS and make the 
averaged matrix element smaller. 

The energy of the FM and AFM fluctuations, tko s f, is 
estimated from a Heisenberg model, where the parame- 
ter J, the exchange integral, is the maximum amplitude 
of the spin wavetj. For ferromagnetic-fluctuations it has 
been proposed that hjj s j = 1/(47V5)E3. The correspond- 
ing property in the present approach is J/4. 

In order to make the calculations of superconductivity 
mediated by spin fluctuations analogous to that mediated 
by phonons, one needs to include a number of representa- 
tive g-vectors. We consider FM fluctuations (q=0) and a 
few AFM fluctuations with q at the zone boundary. For 
fee Co we calculate intermediate q by use of supercells 
with up to 8 atoms. The results from the intermediate q 
vectors show no surprises, but the coupling and J are be- 
tween the values for large and zero q. This indicates that 
no drastic effects from the FS are seen, but it can not be 
excluded that effects from nesting at particular g-vectors 
can be missed by the method of supercell calculations. It 
is as for frozen phonon calculations, where one assumes a 
smooth variation between the values calculated at a few 
g-points. If a special mode-softening is known initially it 
can be studied carefully. 

The strong coupling parametrization of the Eliashberg 
theory for electron-phonon coupling by McMillario con- 
tains the weighted averages of the coupling A and the 
frequency w ;og . 

X = 2 f a 2 F(u) — (4) 



uji g = exp(- / a F(uj)dio) (5) 

X J LO 

With the few number of spin waves we replace the in- 
tegration by a summation over g-points "i" with relative 
weights itfj. The density-of-state contribution a 2 F(w i ) 
from mode i is contained in the value Xi, which is de- 
fined 2Aco i a 2 F(uj i )/uJi for a single mode. 

A = 2J w iK (6) 



uji og = exp(- InuiiWiXi) (7) 

i 

Superconductivity based on electron-phonon interac- 
tion is easily destroyed by magnetic impurities and mag- 
netic fluctuations. The spin-splitting of the potential 
leads to an attraction between electrons of equal spin and 
compete with the attractive interaction between Cooper 
pairs. Effects from impurities are not included, but for 
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the effect from spin fluctuations we use the Daams et al 
modificationE] of the McMillan formula due to A s /, so 
that X e ff = A e)9 /(1 + A s /) will replace X ep . 



knT, 



-1.04(1 + A.//) 



(A e// - M *(l + 0.62A e// ))' 



(8) 



In addition, with /u*=0, the expression becomes formally 
very similar to the expression for Tc by paramagnon 
coupling, see below. 

For superconductivity mediated by spin fluctuations 
the situation is reversed, so that non-magnetic impurities 
and strong electron-phonon coupling reduce Tc. This is 
because the pairing is between electrons of equal spins, 
and scattering with impurities or phonons will reduce the 
spin-splitting of the local potential. For estimating the 
superconducting transition temperature we use a simpli- 
fied version of the weak coupling BCS-like formula for 
paramagnon couphngEita 

k B T c = {hu sf /l.2)exp(-(l + A ep + A° / )/A^ / ) (9) 

The calculated enhancement parameter due to spin fluc- 
tuations is mainly considered to be of p-character, A^, 
although we have not investigated the (k, fc')-dependence 
over the FS. The s-wave fraction of the total A s /, A°^, 
is expected to be small, but it can be more important 
at larger q, when the perturbing potential has stronger 
dipolar character, as in the case of atomic displacements 
in electron-phonon coupling. The s-part should also con- 
tain a contribution from transverse fluctuations!!]. Two 
extreme values of \®j are assumed, either that A|L=A^ 
or that in order to obtain the range of extreme 

values on Tc. The second choice turns out to give a Tc 
of a reasonable amplitude compared to what is measured 
for hep Fe, while the Tc is about half as large if the s- 
part is assumed to be equal to the p-part. It is obvious 
that with the various approximations it is not possible to 
give precise values for Tc, but the values should be taken 
as a possible range of Tc in the different materials. The 
result for hep Fe serves as a calibration; since the Tc and 
its pressure variation agree quite well with experiment, 
there is some hope that the same method can be used as 
a guidance for the other materials, especially for fee Co 
where no experimental data are known. 



III. RESULTS. 

The exchange enhancement is defined as the increase 
of the exchange splitting of the d-band divided with the 
energy of the applied magnetic fieldEj. This definition 
is natural for FM cases and it leads to the usual Stoner 
enhancement for paramagnets. We also define a local ex- 
change enhancement for AFM configurations, when the 
field is applied differently on different sites. The results 
show that these exchange enhancements are largest when 
the material is close to, but on the non-magnetic side of, 



a magnetic (FM or AFM) instability. As the volume is 
decreasing in a magnetic system one will find a critical 
volume where the (local) moments disappear at a certain 
volume. The local exchange enhancement will be infi- 
nite at the critical volume and still large for decreasing 
volumes. But within the region where the volumes are 
larger than the critical volume, the exchange enhance- 
ment drops rather quickly. It is as if the minimum of 
the total energy as function of the moment in a mag- 
netic system is deeper than in a non-magnetic one, and 
the additional fields will not have a large influence on 
the moment. In other words, non-magnetic materials are 
magnetically "softer" than the materials which already 
have a magnetic ground state. A consequence of this is 
that A s / will behave asymmetrically as function of pres- 
sure around the critical volume, so that it will be lower 
in magnetic region than in non-magnetic one. This is at 
least the results from the calculations for several config- 
urations in these pure elements. One may add that we 
have not considered more than the three most common 
structures. The calculations can not exclude that other 
structures may appear at high pressure. 



A. Nickel. 

The calculations for fee-, bec- and hep-structures of Ni 
show that the FM state is very persistent. This is for 
pressures varying from zero up to more than 3 Mbars, or 
when the volume is compressed to 45 (a. it. ) 3 /atom from 
the equilibrium volume of 73 (a.u.) 3 /atom. Only the 
bee structure has a vanishing moment for volumes below 
about 62 (a.w.) 3 /atom, but the total energy is always 
larger than for magnetic fee, about 4 mRy/atom within 
the whole range of volumes. The hep structure is FM 
with a moment ranging from 0.6 to 0.49 /xs/atom in this 
range of volumes. This is very similar to fee, where the 
range of m is from 0.62 to 0.48 /j,s/atom. However, the 
total energy of the hep phase is 3-5 mRy/atom larger 
than in fee. The results suggest that NM bee should be 
stable over FM hep at large pressure, but FM fee is the 
most stable structure everywhere. 

These facts, the stability of the FM fee phase, and the 
low values of exchange enhancements and A s / for mag- 
netic phases, make together that Ni seems not interesting 
as a candidate for superconductivity mediated by spin 
fluctuations. The NM bee phase could have large X s f 
near the magnetic instability for a pressure of the order 
0.5 Mbar, but since it is not the most stable phase we 
refrain from making detailed calculations. 



B. Iron. 

Several results for Fe in a previous publication^ are 
complemented by adding a few GGA results. The 
FM bee phase is stable for volumes larger than 72 
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(a.u.) 3 /atom, while NM hep is stable for lower volumes, 
see fig. 1. Two AFM configurations are stable for hep. 
The AFM-I configuration (when the moments alternate 
along the z-axis) is stable only at large volumes, while 
the AFM-II configuration (alternation along the x-axis) 
is stabilized just below the critical volume where bec is 
stable. This is in agreement with independent calcula- 
tions, where interpretation of iiaman shifts lends support 
to the existence of AFM ordeiH. However, the moment is 
small and the total energy is almost the same for the NM 
hep structure for volumes smaller than 72 (a.w.) 3 /atomtj. 
The calculation of X ep is done in the rigid WS approx- 
imation with the phonon denominator scaled by Debye 
temperature and the P-dependence of the bulk modulus. 
The X ep for iiie-bcc structure agree with previous RMTA 
calculationsEaGa. The values for the hep phase are some- 
what smaller than in ref B, with a similar relative decrease 
as function of pressure (P). This is mainly because of the 
difference in force constant K, while the matrix elements 
for the electron phonon coupling are similar. 

The calculations of X s f and T c are done following the 
procedures described above using the previous results 
from FM and the two AFM configurations^]. The weights 
in the summations for X s f and u>i og (1/5, 1/5 and 3/5 for 
FM, AFM-I and AFM-II, respectively) are obtained from 
the q- values of the spin configurations. The results are 
given in Table 1. The variation of X ep is continuous be- 
tween the largest and smallest volume. This is true also 
for X s f and ui for FM fluctuations. The FM instability 
occurs well beyond the range of volumes considered here 
where the Stoner enhancement never is larger than 3.2. 
The corresponding values A s / (and ui) for the two AFM 
configurations show pronounced peaks (dips) near the 
critical volumes for the magnetic instabilities. The local 
Stoner factors diverge at these volumes. These instabili- 
ties are reflected in the averaged values of X s f and uji g 1 
so that together with the contribution from FM fluctua- 
tions there are conditions for superconductivity near two 
volumes, see fig. 2. The calculated Tc agrees well with 
experiment when X° 8 f is assumed to be zero, with a max- 
imum of about 1.2 K just below the volume at which the 
hep phase becomes stable and with Tc approaching zero 
at larger pressures. A second peak in Tc is seen at larger 
volume, where the AFM-I becomes unstable. However, 
this can not be tested experimentally since the bec struc- 
ture is stable at such volumes. Also hypothetically, it is 
possible that strong FM fluctuations may give supercon- 
ductivity at even larger volumes. As a comparison we 
also calculate the Tc based on electron-phonon coupling 
with reduction due to X s f , using the formula of Daams et 
aM With a**=0, T c is about 1.0 K just after the bec- hep 
transition at the volume 70 (a.u.) 3 /atom. This value is in 
reasonable agreement with experiment, but the reduction 
with pressure is too slow in comparison with experiment. 
At 58 (a. it. ) 3 /atom Tc is still 0.3 K. This corresponds 
to a rather smooth reduction of Tc from 1.0 to 0.3 K 
within a pressure range of ~ 700 kBar, compared to a 
confinement of Tc within 200-300 kBar from spin fluc- 



tuations and within ~ 150 kBar experimentally!!!. Thus, 
from the pressure dependence of Tc it can be concluded 
that superconductivity is more likely to be mediated by 
spin fluctuations. A similar conclusion was reached by 
Mazin et aM 

The good agreement with experiment for the ampli- 
tude of Tc could be accidental. Many approximations 
are used and the statistics for the q-average is poor. 
On the other hand, the relative variations with pressure 
should be more reliable, since the procedure for calculat- 
ing Tc always is the same. This motivates a study of 
Co and Pd with the same procedure in order to verify if 
Co could be superconducting at high pressure and if fee 
Pd, with its large exchange enhancement, will incorrectly 
appear as a superconductor by this theory. 

C. Cobalt. 

The stable phase of Co at zero pressure is FM hep 
with a moment of 1.59 ^a/atom when volume is 73.1 
(a.u.) 3 /atom. The total energy of the NM hep state 
is 15.5 mRy/atom higher at the same volume, which is 
close to the experimental volume at equilibrium. Figure 
3 shows the total energies from FM and NM GGA cal- 
culations of the three structures, and the FM moments 
are shown in fig. 4. The bec structure is magnetic within 
a wide range of pressures, but its total energy is con- 
siderably above those of fee and hep, and the bec phase 
is not of interest here. It is seen in fig. 3 that the fee 
structure becomes stable (over hep) at a volume near 58 
(a.u.) 3 /atom. The pressure at this volume is about 0.8 
and 0.7 Mbar for hep and fee, respectively. The transi- 
tion volume is below the volume (about 60 (a.u.) 3 /atom) 
where fee becomes NM. Therefore, from the experience 
with Fe it is expected that the Stoner enhancement and 
X s f for FM fluctuations become very large at volumes 
just below the FM instability. In spin polarized calcu- 
lations near the critical volume for the FM instability 
it is possible to find metastable low- and high-moment 
states with almost degenerate total energies. This leads 
to large non-linear Stoner enhancements so that S are 
large at large field. The S-values vary from 4 at the vol- 
ume 54 (a.u.) 3 /atom when the applied field is 1 mRy to 
28 at 59 (a.u.) 3 /atom for the field of 2 mRy. As was dis- 
cussed earlier, it is generally easier to find large A s / for 
FM than for AFM fluctuations^, and here we obtain very 
large values, from 0.9 to 2.4 within this range of volume 
and field. This is much larger than for FM fluctuations 
in hep Fe, where the FM instability is outside the consid- 
ered range of volumes. The AFM instabilities are within 
the range volumes, but the corresponding A s / is rather 
modest, since the matrix element for AFM fluctuations 
is not so efficient. 

The FM Stoner factors on the magnetic side of the 
transition are much smaller, about 1.7, and typical A s / 
is of the order 0.1. This agree with the findings for AFM 
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fluctuations in Fe, and FM fluctuations in hep Co. The 
latter is magnetically ordered within the entire range of 
volume where it is the stable phase, and it is not promis- 
ing for superconductivity. 

A few AFM configurations are made for two volumes 
of fee Co just below the volume where FM disappears, 
at 56 and 59 (a. it. ) 3 /atom. The applied field is used to 
impose AFM modulations in the x-direction for different 
wave lengths, corresponding to 2, 4 and 8 atoms in su- 
percell calculations. (The field amplitude, 1 or 2 mRy, 
is the same on each site so that the wave is step-like 
rather than sinus-like.) The local Stoner enhancements 
are lower than for the FM case at equivalent volumes. 
For the shortest AFM wave S is 2.3 compared to 28 for 
the FM configuration. For the longest wave there are 
4 layers of atoms where the field is of the same sign. 
Thus, the two in the middle are surrounded by sites of 
the same polarization and S is 10 within this region of 
"local" FM order. But the two atoms at the edge of 
the 4-layer atoms, have one neighbor with reversed po- 
larization, and S is lower, about 6, for these sites. The 
coupling A s / is relatively small for all of the AFM config- 
urations. It varies from 0.26 to 0.09 and to 0.03 for the 
longest, the intermediate and the shortest wave, when 
the volume is 59 (a.u.) 3 /atom. These values are reduced 
further to 0.07, 0.06 and 0.03, when the volume is 56 
(a.w.) 3 /atom. 

Some values the electron-phonon coupling were calcu- 
lated for the hep phase. At large volume this phase is 
FM and the coupling is dominated by the minority band 
which has the largest DOS. At the equilibrium volume 
X ep is 0.19 for the minority band and 0.05 for the major- 
ity band, while at 54 (a.u.) 3 /atom, when the magnetism 
is strongly reduced, the amplitudes are reversed, 0.25 and 
0.20. This is because of a large increase of the majority 
DOS. A X ep of this order should make superconductivity 
possible if spin fluctuations are absent. However, hep Co 
is still magnetic at this volume and the fact that fee is the 
stable phase at this volume makes it not very interest- 
ing to study superconductivity based on electron-phonon 
coupling. Instead, we calculate A ep for the NM fee phase 
in order to investigate whether it can give superconduc- 
tivity mediated by phonons or if it will prevent supercon- 
ductivity mediated by spin fluctuations. The calculated 
values of X ep are rather stable as function of volume, 
about 0.41. This is for a Qjj of 550 K at the volume 
56 (a.M.) 3 /atom. The calculation of Tc is made as for 
Fe. If the reducing effect from the (q- weighted) A s / is in- 
cluded one obtains 0.1 K at 56 (a.u.) 3 /atom and 0.01 at 
59 (a.u.) 3 /atom. Without influence of spin fluctuations 
we obtain 0.56 and 0.3 K. 

The estimation of Tc* from spin fluctuations are done 
as in the case of Fe, using the q- weighted results in table 
2, with the reducing effect from electron-phonon inter- 
action and with two assumptions for £ — part of A s / 
(either or equal to the £ — 1 part). The results are 
0.2 and 0.06 K, respectively, at the smallest of the two 
volumes, while 1.1 and 0.4 K, respectively, at the large 



volume. 

There are many approximations and uncertainties be- 
hind these values, but since the method is identical to the 
method used for Fe, where the experimental amplitude 
is fairly well reproduced or even underestimated, there is 
some hope that superconductivity could be observed in 
Co just after the hep-to-fee phase transition. The calcu- 
lations predict the transition pressure to be very high, of 
the order 700 kBar, but as seen in fig 3, a small error 
in the relative position between the total energy curve 
will shift the transition volume (and pressure) a lot. If 
the critical pressure in reality were even higher, it would 
be less likely to expect superconductivity mediated by 
spin fluctuations, because the material would be further 
away from the instability point when the fee structure 
becomes NM. The critical volume for the FM to NM 
transition is not as sensitive because if does not depend 
on the relative differences between the total energy of 
two structures. However, the magnetic transition for the 
fee phase is diffuse and extends over some interval of vol- 
umes, where non-linear effects and metastable states are 
found. Therefore, if the hep to fee transition turns out to 
happen at larger volume, there is some margin for having 
larger values of A s / and Tc, than what has been calcu- 
lated here. This is since the volume for the onset of FM 
in the fee structure appears to around 60 (a.u.) 3 /atom, 
while the largest volume for the calculation of A s / is 59 
(a.u.) 3 /atom. All partial values of the enhancements and 
A s /, FM or AFM, tend to increase when the system is 
closer to the critical volume, but non-linear effects make 
the supercell calculations more difficult. 

Some differences with Fe can be noted. The structures 
are different and it is difficult to say if this somehow will 
annihilate the advantages coming from the fact that the 
method of calculation is as identical as possible for the 
two systems. The relative weights between FM and AFM 
fluctuations are also different. Fe is far from the FM 
instability, but there two stable AFM configurations are 
close and fluctuations from them contribute to the final 
X s f . No AFM instability is found in these calculations for 
fee Co. It is the FM case which dominates, with a Stoner 
factor which at the largest volumes becomes considerably 
larger than in Pd. Since Pd at ambient pressure also is 
of fee structure, has a large exchange enhancement, but 
has not been found superconducting, one could test the 
same method for Pd. 



D. Palladium. 

For Pd we consider only the fee phase at the ex- 
perimental lattice constant when the volume is 98.5 
(a.u.) 3 /atom. The clccton-phonon coupling is calculated 
to be 0.36 for a O^ of 274 K. The Stoner enhancement 
for FM configuration iS-4-7, in agreement with earlier es- 
timations using LSDAo. This enhancement is smaller 
than the largest values calculated for Co, which implies 



6 



that the hypothetical FM instability in Pd would require 
a considerable increase in volume. Also the local en- 
hancements for the AFM configurations are smaller than 
the corresponding values for Co. From the calculations 
using 2 and 4 atoms/cell S is 1.7 and 1.9, and for the 
two inequivalent sites in the 8 atom cell the values are 
2.4 and 3.6. The coupling X s f for the three AFM cases 
are of the order 0.02 to 0.14 from the shortest to the 
longest wave, while for the FM case it is calculated to be 
between 1.4 and 1.8 depending on the amplitude of the 
applied field. Thus, the enhancements in Pd is compara- 
ble to the result for Co calculated for the smallest volume. 
The q-averaging is done for the equivalent modes as for 
Co, and the results are A s /=0.15 and oj—13 mRy. 

The estimate for Tc is below 0.01 K for superconduc- 
tivity mediated by electron-phonon coupling (but limited 
by spin fluctuations), and between 0.08 and 0.03 K for 
the two assumptions of in the calculation of Tc from 
ESP (and limited by electron-phonon coupling). The lat- 
ter values are smaller than the results for two volumes of 
fee Co, and Tc is very low despite the large exchange 
enhancement of Pd. The reason is that the large ex- 
change enhancement is not extending to large q- values, 
and the averaged A s / is not large enough for supercon- 
ductivity. The material does not support spin excitations 
where adjacent atoms have different polarization. This is 
very different from hep Fe where opposite polarization on 
near neighbors even give stable AFM configurations. If 
the volume of Pd were ~ 5 percent larger, one could ex- 
pect results similar to what is calculated for fee Co. The 
reduction of enhancements at large q is crucial for an 
understanding of the small Tc or even absence of super- 
conductivity. A calculation based only on the standard 
Stoner enhancement (at q=0) would not be correct. A re- 
cent work suggests that_spin-orbit coupling can suppress 
superconductivity in Pto. One may speculate that such 
an effect should be more evident in Pd than in Co because 
of the relative strengths of the spin-orbit coupling. 

IV. CONCLUSION. 

Two possibilities for superconductivity in four ex- 
change enhanced metals have been investigated by cal- 
culations of the coupling constants for electron-phonon 
coupling and spin fluctuations. The method for calcu- 
lating \ ep and the corresponding Tc has been quite suc- 
cessful in many studies during the last decades, while 
the method for calculating A s / is new and needed some 
testing. Both methods include approximations and one 
should not rely too much on to the absolute values of 
Tc, but look more at the relative variations among the 
different cases and different pressures. Despite these 
precautions, it seems clear that superconductivity based 
electron-phonon coupling is unlikely in Pd, Co and Ni. 
This is because of the limitation of Tc coming from spin 
fluctuations or from stable magnetic moments, and from 



the rather modest values of X ep . Hep Fe could have a rea- 
sonable Tc, but the agreement with experiment is poor 
for the pressure dependence when electron-phonon cou- 
pling determines Tc- 

The results for Tc based on ESP are partly a test 
of the method, partly a prediction. The experimental 
situation in hep iron is surprisingly well reproduced, both 
the amplitude of Tc (~ 1.2 K) and its variation with 
pressure. This result indicate that a representative q- 
point average of A s / has to be calculated. The calculation 
for fee Pd pass the test as well, since Tc is very low (~ 
0.08 K) despite its large Stoner factor. Again, this is 
the result only when several configurations are weighted 
together. The calculations indicate that a NM fee phase 
of Co becomes stable at a high pressure of the order half 
a Mbar. Here, Co behaves as an improved version of 
Pd with larger enhancements, but with similar relative 
variation as function of q. The improved conditions for 
superconductivity from spin fluctuations give a Tc of the 
order 1 K. 

These calculations do not include considerations of de- 
fects. Since superconductivity due-, to FM fluctuations 
is sensitive to impurity scattering^, it is possible that 
lattice defects, induced in pressure experiments near a 
structural transition, can suppress Tc- It is not yet clear 
if experimental conditions can be improved to make Tc 
higher in hep Fe or even in Pd. The reduction of Tc 
compared to Tc itself is proportional to vp/(£A), where 
I the impurity scattering length, A the superconducting 
energy gap and vp is the Fermi velocityEj. Knowledge 
of the q-dependence of the destructive effect from defects 
on A s / or ESP should be useful, since this can make a 
difference between Fe and Co (and Pd). The FM fluctu- 
ations are relatively important in the latter two metals 
of fee structure. This is in contrast to hep Fe, where 
the AFM fluctuations have a large weight. Calculations 
for supercells containing a defect give information about 
how much a defect will hamper the development of FM 
or AFM spin waves. The real space interpretation of the 
AFM results for Pd is that it is not easy for neighboring 
sites to have different magnetization. If so, it should also 
be difficult for a FM moment to develop near to a non- 
magnetic impurity in Pd, which could make this material 
more sensitive to impurities than Fe. In the absence of 
full calculation with impurities one can also note that 
from the dispersion of spin waves it follows that fluctu- 
ations with intermediate q are mobile with large values 
of the group velocities. They may interfere with defects, 
while waves with dut/dk ~ (as for a purely FM fluc- 
tuation) can develop easier between defects. Co is an 
intermediate case between Fe and Pd even without con- 
siderations of defects, still with possibilities for a mea- 
surable Tc- 
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TABLE I. Volume per atom, partial and q-weighted X s f 
and u) s f (in energy units) for hep iron. The weights are 0.2 
0.2 and 0.6 for FM, AFM-I and AFM-II, respectively. The 
electron-phonon coupling decreases continuously from 0.37 for 
the largest volume to 0.14 for the smallest volume. 
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3.5 
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0.15 


0.02 


0.06 


0.07 


11. 


6.2 


3.2 


5.7 



TABLE II. Volume per atom, A ep , partial and total 
q-weighted A s /, and total u>i og (in energy units) for spin fluc- 
tuations in fee Co at two volumes, and in fee Pd (last line). 
The weights are 1/16, 3/16, 6/16 and 6/16 for FM, and the 3 
different AFM configurations ordered according to the wave 
length (2 4 or 8 atoms). 
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FIG. 1. Calculated total energies for paramagnetic hep 
(heavy line) and ferromagnetic bec (thin line) iron using the 
GGA potential. 
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FIG. 3. Calculated total energies using GGA potential as 
function of the volume for hep, fee and bee Co. 
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FIG. 2. Calculated Tc in hep iron from eq. |^ using the 
values A ep and the weighted values of A s / and u>i og in Table 1. 
The full line is for = and the broken line for \°f — Ajj. 
To the right of the vertical line is the stable structure bee (and 
Tc=0) according to the calculations. Superconductivity has 
been observed with Tc below 2 K in a pressure range oij-about 
150 kBar starting near the bec-hep transition pressureEJ. 
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FIG. 4. Magnetic moments calculated using GGA poten- 
tial for hep, fee and bee Co. 
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